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1 Introduction 



The low dimensional {d = 0, 1, 2) integrable models, for a long time a sub- 
ject of interest in the mathematical physics, are now attracting much attention 
also in connection with problems in condensed matter physics. The develop- 
ment of nanostructure technology nowadays allows us to consider the problems 
in zero(quantum dots), one(quantum wires, organic polymers) and two(hall effect, 
high temperature superconductivity) spatial dimensions. 

The well known one dimensional Heisenberg model of a nearest-neighbour 
interacting chain of | spins, 

TV N 

H = J2H,,+, = J2a,a,+,, N + 1 = 1, (1) 

i=l 1=1 

was solved in 1931 by Bethe |^ using what is now called the Bethe Ansatz. 
During the last 25 years this method has been developed into what is now called 
the Quantum Inverse Scattering Method(QISM) 0, ^ (or the Algebraic Bethe 
Ansatz(ABA)). The method allows one to find a complete set of eigenvalues and 
eigenstates of the Hamiltonian of a one dimensional integrable model. 

The origin of integrability is a set of equations, first derived in the article by 
Yang [H, and which, after the works of Baxter 0, appeared to be equivalent 
to existence of an infinite number conservation laws in the model, and called 
Yang-Baxter Equations (YBE) as a consequence. 

However, the ABA technique was mainly developed in the context of the XX Z 
spin chain model and its application to models formulated originally in fermionic 
operator language was complicated. For example, in the Hubbard model, solved 
in the very beginning by a Coordinate Bethe Ansatz (CBA) method 0, the 
implementation of ABA to find a full solution was made only recently 0] after 
remarkable step done by Shastry [0! 01- He used the Jordan- Wigner trans- 
formation to pass from fermionic to spin language, in order to apply QISM, 
formulated in matrix form. 

Recently the interest of formulating spin chain models in terms of Fermi op- 
erators by use of Jordan- Wigner transformation []TU|, [TI|, and by use some new 
method developed in ||T2|, Il3|, 113, has increased. It appeared that the explicit 



construction of integrable models and application of QISM in terms of Fermi op- 
erators becomes easy. Because the Fock space of fermions is finite dimensional the 
investigation of the complete set of eigenstates and eigenvalues in this language 
becomes relatively simple. 

It is well known that the trivial generalisation of the Heisenberg model to 
higher spins is nonintegrable and gap-full |T^. The first integrable model for the 



spin-S* was established by G.Uimin |15] and later by J.K.Lai |16] and B.Sutherland 
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They have considered the Hamiltonian 



N 

i/ = J^(M,,+i + P,,,+i), N + l = l, (2) 

1=1 

where is the operator permuting spins at the sites i and i+1 and Idi^i^i is the 
identity operator, shown its integrabihty and found a solution by the Coordinate 
Bethe Ansatz. Another integrable model for arbitrary spin-S* was constructed in 
the articles [TO, 120]. As it appeared, the integrable Hamiltonian is a polynomial 



Q2S of degree 2S of the product of nearest-neighbour spins S'jS'j+i 



TV 



Hs = JJ2Q2s{S^S^+l) N + l = l, (3) 

2=1 

where 

Q2.(x) = -E(E)nfTV' 

j=l k=l 1=0 -^j 

a;, = ^[/(/ + l)- 25(5+1)] (4) 

and has a critical point, near which the model is equivalent to Wess-Zimino- 
Witten-Novikov (WZWN) model with corresponding coupling constant k = 2S 
2^ . However in the article ]jl9|] only 2^ (L-is the length of the chain) states from 



all 3^ amount of states were found together with corresponding Bethe equations 
for them. 

The full set of SU(2) invariant, integrable, isotropic quantum spin chain mod- 
els for spins S <6 was found in |2T[] . 

In this article we will use the technique developed in [l^, which is alter- 



native to the Jordan- Wigner transformation, to fermionize the integrable spin-1 
ULS isotropic chain model and find its Hamiltonian in terms of fermionic fields, 
which will demonstrate its equivalence to integrable supersymmetric (SUSY) t—J 



model ]]23[. After that we generalise this construction to spin-S" and carry out 
the fermionization of the spin-S" ULS chain model. The resulting Hamiltonian 
can be considered as spin-S* analogy for the t — J model. 

The article is organised as follows. In the Section 2 we describe briefly the 
fermionization scheme for integrable models and QISM. 

The fermionization of the spin-1 Uimin-Lai-Sutherland (ULS) model will be 
carried out in the Section 3. The Hamiltonian will be written in Fermi terms and 
coincidence with SUSY t — J model will be shown. 

In the section 4 we will construct the fermionic representation of the SU{2) 
algebra for spins (25* + 1) < 8 and apply that construction to find the i?-matrix 
and corresponding Hamiltonian in terms of fermions in section 5. 
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2 Yang-Baxter equation and structure of the 
quantum spaces of the integrable model. 

The key of integrability of the model is the Yang-Baxter equation (YBE), which 
imphes some restrictions on Raj-Taathx - the constituent of the transfer matrix 
T{u), defined as 



N 



(5) 



T (u) ^ trail Raj (U) . 

The YBE ensures the necessary and sufficient conditions of commutativity of the 
transfer-matrices for different values of the spectral parameter u, which has a 
physical meaning of the rapidity of pseudoparticles of the model 



[T{u),T{v)]=0. 



(6) 



The Hamiltonian of the model can be obtained by logarithmic derivative of the 
transfer matrix at zero spectral parameter 



H 



d\nT{u) 
du 



\u=0 



(7) 



By definition Raj acts as a intertwining operator on the space of direct product 
of the so called auxiliary T4 (v) and quantum Vj (u) spaces 



Raj (U, V) : Vaiu) ® Vj{v) Vj{v) ® Vaiu) 

and can be represented graphically as 



(8) 



Raj(u,v) =■ 



Va{u) 



Vj{v) 



Fig.l. 



The spaces Va{u) and Vj{v) with spectral parameters u and v are irreducible 
representations of the affine quantum group UqCj, which is the symmetry group 
of the integrable model under consideration. Provided that the states | a) e \4 
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and I j) & Vj form a basis for the spaces Va and Vj, following [T^, we can 
represent the action of the operator Raj as follows 



Raj \3)®\a)= I a')® I f), 



(9) 



where the summation is over the repeating indices a' and j' (but not over a and 
By introducing the operators 



a') {a I; 



f){j 



in the graded spaces Va and l^- correspondingly, one can easily rewrite 
follows 



(10) 
1 as 



Ra, = Ra, I j) I a){a I (j 1= I «') I I ij 1= (-l)^^'^^^^^"^ K/^^'^i' 

(11) 

where the sign factor takes into account the possible grading of the states | a) 
and I j), p{a) and p{j) denote the corresponding parities and the summation over 
the repeating indices. It is easy to see from the expression (0) that the operators 
Xa have the projection property 



In the conventional form the YBE looks like 



i?^?,, {u, v)R% {u, w)B^^.„ {v, w) = B^^, {v, w)R%„ (n, w)Ri!!l„ (n, v) 



(12) 
(13) 



By use of equations (0) and (|T^ the matrix-valued YBE eq. (0) one can easily 
be transformed to the following operator valued equation 



Rahiu, v)Raj{u, w)Rbj{v, w) = Rbj{v, w)Raj{u, w)Rab{u, v) . 



(14) 



Using the graphical representation of the -R-operator presented in Fig.l, one 
can draw the YB equation graphically as follows 



1 




Vj{w 


Va{u) \^ 


Va{u) 




Vaiv) 


Vaiv) 






Fig.2 



An alternative form of the YBE above is the commutativity of the i?-operator 
(|TTp with the co-product A of the corresponding quantum group Ugg |T^, |TB . 
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Let us introduce now the L-operator by taking matrix element of the R- 
operator (|ll]) in the auxihary space Va with the basic vectors | a) 



L:, (u, v) = {a' I Riu, v) I a) = {-if^^^^^^ E^^, (u, v) Xj' (15) 
Then, following Jl4|, it is easy to show that YBE takes the form 

Rtr (L ®s irj = {L ®s L)i'' (16) 

where -RJJ^ = R^'^ and graded direct product is defined as follows 

{A ®, B)Z = (_i)W'^)+fW)fW A^B^a, (17) 



which is conventional in the graded Quantum Inverse Scattering Method |25 



3 The spin 1 ULS-chain in terms of fermions 

Now we would like to apply the approach, described above, to fermionize the 
spin 1 integrable isotropic ULS chain with the Hamiltonian (|^). 

To begin with we should realize the spin 1 algebra in terms of fermions. At 
least two sort of fermions are needed in order to express three basic states | +), 
I 0), I — ) of the spinl particle with the z component of the spin equal to 1, 0, — 1 
respectively. Let us define c^, c^, cr =t, | as a creation -annihilation operators 
of fermions with the up and down spins correspondingly, together with their Fock 
space I 0), I 0"). 

The generators of spinl representation of the 5^(2) algebra can be built by 
use of this Fermi operators as follows: 

S+ = (l-nt)q + (l-ni)c+, 

S- = (l-nT)c+ + (l-nJct, (18) 
S, = n^-ui, 

where Ua is the fermion number operators. 

By this definition of su{2) generators, the states with definite third projection 
are realized through fermionic Fock space as follows 

|+)=|T)|0), |0)^|0)|0), |-)=|0)|i). (19) 

Indeed, easy to check that Sz\+) = |+), •S'^lO) = 0, Sz\—) = — |— ), 'S'_|_|+) = 0, 
5+|0) = |+), = |0), = |0), 5_|0) = S^\-) = 0. 

As it is obvious from formulas ([T9|), we have constructed a graded space with 
the following parities for the basic vectors 

Pi\ +)) = p{\ -)) = I, Pi\0)) = 0. (20) 
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The sum of the square of spin generators, which is the Casimir operator 

S+S- + S-S+ + Sl=C = 2{l-n^n^), SX = 2Sa = CSa (21) 

determines the projection operator 

U=l-n^ni, {1 - n^ni)^ = 1 - n^ni, (22) 

which maps 4-dimensional space of direct product of two spinl/2 spaces onto 
3-dimensional subspace spinl. |C acts as an identity operator on eigenstates of 
the 5*2: 

-C \ m) =\ m) m = +,-,0, (23) 



and annihilates the fourth state 



C|T) |i) = o, 



(24) 



which corresponds to scalar state in the expansion of direct product of two half 
spins. 

In the spinl case the operator X takes the form 





1 


|-)(-| 




|-)(o| 


|-)(+l \ 


\m)i{k\i = 




|o)(-| 




|0)(0| 


|0)(+| 




\ 






l+)(0| 






1 


(1-^1 




(1 - n^)c+ 








(1-^1 




(i-^t)(i- 


ni) C|(l-n|; 




\ 






4(1 -nj 


^t(i - 



(25) 



where i denotes the chain site. Trace of this operator is an identity operator due 
to the completeness of set of the states. 

The operator Uij, which permutes the states between spaces Vi and Vj has 
the form 



rrii 



(26) 



The sign in this expression reflects the grading of the spaces Vi 

One can find the i?-matrix of the spinl model in its bosonic form in |1^, |T6|, ^ 
but it can easily be guessed also from the expression of the Hamiltonian (0), 



^ i=l 



(27) 



which can be represented also as Unity + Permutation operators in the 3-dimensional 
spinl space. It has a following form 



U 



6^51 + u 5^151,, m, k,s,q = +, 0, 



(2^ 
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where first term is unit and second term bosonic permutation matrices. 

The L-operator of the model, defined in fermionic representation by the for- 
mula (|15|), takes the form 



LJu) 



f M - (1 - n|) (1 - n{) (1 - n|) 



(1 -^T)ci 
( u + \{Sl- 



2 

SzS^ 



—S^Sz 



u 



\ 



Si 

s s 

u+\{Sl + Sz) ) 



(29) 



and the Yang-Baxter equation (|TB|) takes form 
' (v)L'l(u) + (u 



V] 



s ku) Ll^ {v) 

^Y(q)p(s)+pim)p{q)+p{s)p{m) j^q ^^^^ jk 



(30) 



lYioui the equations above one can be easily read off the invariance of the 
YBE under the transformations 



(31) 



where M and are arbitrary numerical matrices of the same grading as L (u). 

According to the prescription of the previous section, by use of formulas (|TTD , 
( PSD and the definition (|^ as a logarithmic derivative of transfer matrix, it is 
easy to find the following expression for the Hamiltonian via graded permutation 
and identity operators 



(32) 



N N 
i=l i=l 

due to the fact that Rik (0) = Pik- 

In terms of Fermi operators this Hamiltonian can be represented in very simple 
way as 

= AiAj ihj + Vij,^V^j,i) Ai A, , (33) 
where Aj is projection operator defined in (p2D and Vij^a = 1 — {cfcr — cf^cr){(^i,(T " 



0" =t, I are the fermionic permutation operators. 



It is now straightforward to recognise in the expression (p3D the Hamiltonian 



of supersymmetric t — J |13, 14, B3, 24] model 



The symmetry (|3T| ) mentioned above is generated by the operator 



( (1 -n|)ra| (1 -^t) 



ib 



(l-^T)ci (1 -n|) (1 -nj,) C|(l-n|) 



Q 



N 



1 (i-^i) 

[H, Qt] = 0. 



(34) 



1=1 
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In fact Qi coincides with the operator Xj (|25|). Being defined as global quantity 
along the chain, the matrix Q is inert under permutation of the site indices, i.e. 
it commute with the local Hamiltonians -f/jj+i. 

In the unrestricted Fock space of spin-up and spin-down fermions one can 
construct the X-operator (p!OD, which will be a graded direct product of two X's 
in the one-fermion space 

{X)i,'' = {X^®sX{)i,'\ a, 6 = 0,1. (35) 

It is easy to see, that the spin-1 X-operator, defined as in eq.(^), can be obtained 
from the expression (^) by deleting the row and column, corresponding to the 
exclusion of state If) ||). 



4 Spin-S Representation of the SU(2) algebra 
in the Fock space of Fermi fields 

In order to describe the 25+1 dimensional space of the spin — S representation 
by fermions we will consider the Fock space of r copies of the fermionic creation- 
annihilation operators c^, c^, = 1, 2...r. 

V. = (g)V;, (36) 

where is the Fock space of /i-type fermions. The dimension of this space is 2*" 
and hence the necessary condition to have a enough degrees of freedom for the 
spin — S states is {2S + 1) < 2^. 

As usual let us mark the basic elements of the fermionic Fock space by their 
filling numbers = 0, 1; /i = 1, ...r 

Im, ...nr) = |ni) ■ 1^2) ■ ■ ■ \nr) = (c+)"H4)"' " " " (c+)"10) (37) 
\ni,...nj.) E Vr 

It is necessary to order this basic vectors by integers which can be achieved, 
for example, by use of binary system of numbers. But let us define m as a map 
of the set of sequences (ni,n2, ...n^) onto integers 1,2, ..2^ 

m : (ni,n2,...n^) ^ 1,2...2'' (38) 

Then it is clear, that 

Vr=^Wm (39) 

m=l 

where Wm is the space defined by the basic vector \ni,n2, ...rir). 
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There is a natural grading in the Fock space of fermions. Zero and one particle 
states in define a graded homogeneous subspaces = Vo^©Vi^ with the parity 
p as a function Vjfj, — * Z2 

Pilrif,)) = n^, \n^) G K^,;,, n^ = 0,l (40) 

Grading of the tensor product space Vr = <S>J^=iV'^ is defined by parities of its 
constituents, namely, by parities of the basic elements of V,-, as follows 

r 

p{\n,,...nr)) = Y.PiK))- (41) 
Now, following section 2 let us define the operators 

(pK)+pK))Ep(";) 

^™ = X:i:::::; = \n,,...nr){K,...n[\ = i-l) •• 

. . .(^_]_^(P("2)+P(n2))p(n'i)j!^n; _ _ _ X™^, (42) 

where 

xi = K){n'^l /i = l,...r (43) 

is the corresponding X-operator in the Fock space of the /i-type fermions. 

n' 

2^ can be considered as a matrix 



n 

Operators Xn^ can be considered as a matrix X^ with the operator valued entities 



of the form 



= ( ? 1 ■ (44) 



By use of definition ( |17D the expression (|42|) can be represented as a graded 

n' 

direct product of the operators X„^ 

^S::;^ = (^i®.---®.^.):i;:£ (45) 

It is now meaningful to define also the operators which projects an arbitrary 
state from the representation space Vr on to the particular one-dimensional sub- 
space defined by the state |ni, ...rir). Obviously it is 

Pm{ni,...nr) Pni,...nr |Til , . . .Tlr) (^r ; • • •'^l | 

Pm{ni,...nr) ~ Pm{ni,...nr) (46) 

Therefore the operator 

^m{ni,...nr) ^ Pn\,...nr (47) 

excludes the one dimensional subspace, defined by the basic vector ...n^), 
from the whole Fock space Vr, reducing the dimension by one. Hence, in order 
to construct a (25* + 1) dimensional representation space V25+1 for the spin — S 
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we should act by corresponding amount of projectors of the type ( ^71 ) on the 
whole 2''-dimensional space Vr 

V25+I = Ai ■ ■ ■ A2r_(2S'+l)Vr. (48) 

The definition ( ^Sf ) means that any operator A which acts on the space Vr, 
will have a form 

As = Ai ■■ ■ A2r_(25+i)^Ai ■ ■ ■ A2'-_(2S'+i) (49) 

on the space V(2s+i). 

Again, let's now numerate the remaining basic directions in V(2S' + 1) by 
/(rii, ...rir) = 1, 2. ..(25' + 1) and denote the corresponding X-operator as 

x/ = Ai ■ ■ ■ A2r_(25+i)-^^i\':::";Ai • ■ • A2r_(25+i)- (so) 



Then, as it is easy to check by use of formula (p!2|) and some simple algebraic 
calculations, that the operators ,B = 1, 2, 3 defined by 

5^ = (5^)'-X/„ /,/' = !,. ..(25 + 1), (51) 

fulfill the commutation relations of the S'f/(2)-algebra generators. In (|5lD (^S^^^ 
are the ordinary number valued matrix elements of the spin — S generators of 
the SU (2)-algebra, the only nonzero elements of which are 

{S^t' = {S^)l'' = ^{S + M){S-M + 1), 5- = (5+)^ 

(-53)! = {S3)m = M, M=-S + l-l. (52) 

^From now on we will consider the case r = 3 and the spins S < (2'^ — l)/2 = 
7/2 for simplicity. 

For the convenience let us numerate the even elements of basic vectors \ni, ...rir) 
as a, (3 = 1, 2, 3, 4 in the following way 

10,1,1)^1, 11,0,1)^2, |1,1,0)=3, 10,0,0)^4. (53) 
and the odd elements a,b = 5, 6, 7, 8 as 

11,0,0)^5, 10,1,0)^6, |0,0,1) = 7, 11,1,1)^8, (54) 
Hence we have 

4 8 

V3 = W^a = Km © K.en (55) 

a=l a=5 

Now by use of definition (^) and (0) one can find the following expression 
for the -operator (m, m' = 1, ...8) in the unrestricted space V3 
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where 



ctc2n3 



-ctn2C3 



-C1C2 ns 
ni(l - 722)723 



\-{l-ni)c2C3 -Ci(l- 773)03 



C1772C3 
-771C2C3 

77i772(l - 
-CiC2(l 



^1 - ni)cjcj 
ct(l - n2)ct 



-1 

etc 



1 ^2 



-773) (l-77i)(l- 772)(l- 773) / 



,(57) 



/ 77i(l - r72)(l - 773) 

-CiC^(l - 773) 
-ci(l - 772)4 



-402(1-773) 4(1-^^2)03 -^2lC2C3 ^ 

(1-771)772(1-/73) (l-r7i)4c3 C1772C3 

-(l-r7i)c24 (1 - r7i)(l - r72)r73 -0102/73 

-4^24 44% r7i772773, / 



,(5^ 



and 



-4C2C3 -77i(l - 772)03 -77102(1-773) o|(l - 772)(1 - 773) 

-(1-771)77203 010^03 01772(1-773) (1-771)0^(1-773) 

(l-77i)02773 01(1-772)773 -01020;^ (1 -77i)(l -772)0^ 

V 4'^2% -77102^773 771^203" 



444, 



,(59) 



(60) 

It is easy to see from the expressions ( ^T]) and (|58D that X^' can be obtained 
from the by particle-hole transformation of all fermions in a following way 



Cl 



-1 , 



C2 



-2 , 



C3 



(61) 



5 Fermionic representation of the spin-S ULS 
model: Spin-S generalisation of the super- 
symmetric t-J model 

In this section we will fermionize the Hamiltonian of the ULS model defined 
by the formula (||) for spins S < which can be described by three Fermi fields. 
The procedure is based on the approach described in the section 2 and is an 
alternative to the usual Jordan- Wigner transformation. First we consider spin |. 

We start with the i?-matrix of the ULS model 

R'M = + u 777, fc, g = 1, 2, ...8. (62) 

and in the same way as for spin-1 case described in the section 3, we obtain the 
Hamiltonian 

N N N / 3 \ 

Hi =11 H,,i+i = E (^M+i + n^+i) = E + n (63) 

i=l i=X i=l y /i=l J 
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Our goal is now to rewrite this expression of the Hamiltonian in a way which 
is similar to spinl t — J model. 

According to formula (|49|) the Hamiltonian for spins less that | will be 



2 

= Ai ■ ■ ■ A8-(25+l)-f^| A8_(2S+1) ■ ■ ■ Ai, (64) 

where the projector Ai ■ ■ ■ A8_(2S'+i) cuts off some subspace from the full 8- 
dimensional Fock space of representation V3, reducing the degrees of freedom 
to (25" + 1). The different choices of the cutted subspaces (projectors A) corre- 
sponds to isomorphic representations in a mathematical sense and simply means 
particle-hole transformations for some of three fermions from a physical point of 
view. 

Let us now expand the expression ( |63|) for the Hamiltonian with the 

signs, as it should be according to eq.(|ll|) and substitute the expressions for the 
X-operators (|5^ - |60D . Then after some algebra and by use of definitions (|46|) of 
the projectors we obtain 

^.+1 = H]iZ + H^i:^ + H^iZ + Hf^. + ^,+1 - Hf,^r (65) 

The different terms in this Hamiltonian is defined as follows. The terms in H^'^!^ 
are the one and three particle hopping terms 



TTl,hop _ f V \a f Y fV'\"^fV ^" 



= CjA (-Pj,4-Pj+1,4 + PjfiPj+lfi + PjjPj+lJ + Pj,lPj+l,l) 

+ '^t,2 {Pj,5Pj+l,5 + PjaPj+IA + PjjPj+l,! + Pj,2Pj+l,2) Cj+1,2 

+ (-^i:5-Pj+l,5 + PjfiPj+lfi + Pj^Pj+^A + Pj,3Pj+l,3) Cj+1,3 

~ (cj^l9,2Ci,3Cj+l,lC^_l_1^2Cj'+l,3 ~ Cj,iC^2Cj,3c/+i^iCj+i,2C^+i_3 

- Cj,iCj,2C^3C^+i,iC^+i,2Cj+l,3 - Cj^lCj^2C^3Cj+l, 19+1,2^^+1,3 " h.C.]{QQ) 



separate and Hf^?^ hopping terms, correspondingly. 



We have extracted the terms, which contains a = 4 and m = 8 indices into 



H}^ = {x,rAx,^i)t + {x,)Ux,^,)2 (67) 



— ~Pj,S (Cj,2Cj,3Cj-+i + '^i.l'^i.sCj+i^iCj^^i 3 + Cj,iCj,2Cj-^i,iCj+i,2 j Pj+i,s 

As it is clear from the formulas (|59| ) and (|60|) the Hj'^°^ term can be obtained 
from the Hj'^^^ by particle hole transformation ( |6lD These terms are responsible 
for the hopping of the Fermi pairs. 

jjdiag jg ^YiQ diagonal term containing coulomb interaction of fermions 

H''^' = ixxix,+i): + ix,)i{x,^,r^~{x,)i{x,^,): + {j^ 

= (1 ^ Pj,l — Pj,2 ~ Pj,3) {Pj+1,1 + -Pj+1,2 + -Pi+1,3 + Pj+1,4) 

- P,.8(P,.+i_5 + P,.+i^6 + P,.+1^7 + P,.+13) + (j.-^j + l). (68) 
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The last two terms in ( |65D are analogous to the spin-spin interaction term in 




where 



Sf 



Cail - ni){l - n2){l - ns), 



a = 1,2,3. 



(70) 



and {S^)f is the spin — 1 representation of the SU(2)-algebra. 
This expression can be obtained easily from 



(71) 



which follows from formula (pT]). By appropriate an choice of projectors A^, cr = 
1, ...8 — (25*+ 1) we can simplify the expression for the Hamiltonian, for example 
by deleting the terms H^'''"p and H^'''"p. 
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